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We study the dynamical properties of a diffusing lamb captured by a diffusing lion on the complex 
networks with various sizes of N. We find that the life time (T) of a lamb scales as (T) ~ N and 
the survival probability S(N — * oo,t) becomes finite on scale-free networks with degree exponent 
7 > 3. However, S(N, i) for 7 < 3 has a long-living tail on tree-structured scale-free networks 
and decays exponentially on looped scale-free networks. It suggests that the second moment of 
degree distribution (fc 2 ) is the relevant factor for the dynamical properties in diffusive capture 
process. We numerically find that the normalized number of capture events at a node with degree 
k, n(k), decreases as n(k) ~ k~ a . When 7 < 3, n(k) still increases anomalously for k « k max . We 
analytically show that n(k) satisfies the relation n(k) ~ k 2 P{k) and the total number of capture 
events Ntot is proportional to (fc 2 ), which causes the 7 dependent behavior of S(N,t) and (T). 

PACS numbers: 05.40.Fb,89.75.Hc,89.75.Fb 



I. INTRODUCTION 

The diffusion-controlled reactions, in which diffusing 
particles are immediately converted to a product if a pair 
of them meets together, have many physical applications. 
Examples include electron trapping and recombination, 
exciton fusion, wetting, melting, and commensurate- 
incommensurate transitions 0, 0, y, El IE IE ■ Among 
these examples, dynamical properties of wetting, melt- 
ing, and commensurate-incommensurate transition are 
known to be related to the diffusive capture process 
[E El IE IE 0- whose kinetics can be simplified by dif- 
fusing preys-predators model [E III • 

In general, the dynamical properties of diffusion- 
controlled reactions are closely related to the first-passage 
phenomena of random walks (RWs), because the reac- 
tion occurs when diffusing particles first meet a trap or 
an absorbing boundary 0, IE 13 ■ Especially, the diffusive 
capture process can be mapped into coupled RWs , in 
which n predator random walkers (lions or moving trap) 
stalk a prey random walkers (a lamb or a moving par- 
ticle). Since the first-passage probability can be easily 
obtained from the survival probability Q, the survival 
probability of a diffusing particle or a lamb is important 
to understand long time characteristics of the diffusive 
capture process [Ell ■ The survival probability of a lamb 
has been well investigated in d-dimensional hypercubic 
lattice. The diffusive capture process in the spatial di- 
mension d > 2 is unsuccessful and the survival probabil- 
ity is finite for any n. In d = 2, the capture is success- 
ful but not efficient, i.e., the life time of a prey random 
walker is infinite and the survival probability shows a 
logarithmic decay E2- The predators efficiently capture 
a prey for spatial dimension d < 2. On a (i-dimensional 
regular lattice, the survival probability S(t) depends on 
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where the exponent (3 varies as n and the ratio of diffusion 
constant of walks. 

There have been increasing interests on complex net- 
works, since many structures of physically interacting 
systems are shown to form nontriviai complex networks 
El EE EE El EE E3 The early studies on complex 
networks mainly focused on their topological properties. 
Recently, the physical systems whose elements interact 
along the links in complex networks and the effects of 
the underlying topologies on the dynamical properties of 
the system have drawn much attention. Many dynamical 
systems on complex networks show rich behaviors which 
are far from the mean-field expectations pol l2ll I22I |23| . 
Recently, RWs on complex networks have been studied 
to show that the dynamical properties of RWs are closely 
related to the topology of underlying networks |2J, [25| . 
This implies that the first-passage property of the diffu- 
sive capture process is also affected by the topology of un- 
derlying network. Moreover, the several studies on RWs 
have provided very interestin g applic ation to searching 
information over the Internet (26L l27| . For example, the 
major searching engines such as Google and YAHOO use 
general random walking robots along the links between 
hyper-texts to collect information of each web page psj . 
Some peer-to-peer (P2P) network protocols also exploit 
the dynamical properties of RWs on complex networks 
[lEl ■ Especially, the searching algorithm for P2P network 
can be mapped to the system of a diffusing particle to find 
an immobile absorbing particle (static trap). Our recent 
study (23 has also found that the searching algorithm in- 
spired by the diffusive capture process can dramatically 
reduce the traffic over the Internet with scalable search- 
ing time. The diffusive capture process can be also ap- 
plied to investigating condensed phase of conserved-mass 
aggregation model on scale-free networks (SFNs) with 
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2 < 7 < 3 [2jJ . Therefore, it is theoretically and techno- 
logically very interesting to study the effect of nontrivial 
complex networks on the diffusive capture process. 

In order to understand the relationship between the 
non-trivial topology of underlying network and the diffu- 
sive capture process on it, we study the random walking 
lamb and lions on SFNs. The SFNs are the networks 
whose degree distribution follows P(k) ~ fc~ 7 . The de- 
gree exponent 7 of the most natural and technological 
networks is in the range 2 < 7 < 3 0, 1^3] ■ Moreover, 
the probability that a random walker visits a node has 
been shown to be proportional to the degree k of the 
node [24|. So, the probability that two random walking 
particles meet at a node is proportional to k 2 . Therefore, 
we expect that the dynamical properties of the diffusive 
capture process on the SFNs with the degree exponent 
2 < 7 < 3 are different from those on the SFNs with 
7 > 3, because (A; 2 ) diverges for 7 < 3. To verify our ex- 
pectations, we study the diffusive capture process on the 
SFNs. More recently, Noh et al. [25| also showed that 
the return probability of a random walker on the looped 
SFNs (LSFNs) is completely different from that on the 
tree-structured scale-free networks (TSFNs). Therefore, 
the diffusive capture process on TSFNs is also investi- 
gated. 

This paper is organized as follows: In Sec. II, we 
present a method of generating TSFNs and dynamical 
properties of the diffusive capture process on TSFNs. 
The results on LSFNs are shown in Sec. III. In Sec. 
IV, we discuss the relation between degrees of SFNs and 
capture events. Finally we summarize our results in Sec. 
V. 



II. DIFFUSIVE CAPTURE PROCESS ON 
TSFNS 

All the models are defined on a graph Gr = {N,K}, 
where N is the number of nodes in the networks and K 
is the number of degrees with the average degree (k) — 
2K/N. Networks with the sizes N = 10 3 ~ I0 6 are 
generated for each 7. 

To make tree structure, the total number of links 
in networks with the size N is constrained to be 
Kt = 2(N — I). So we use the following method 
to make TSFNs. Select N integer random numbers 
{fci, k2, • ■ • , fc/v} whose probability distribution satisfies 
f(k) ~ fc-T. Let K Q = Eti h- VK T > Kq (K t < K ), 
then choose kj randomly from the set {ki} and decrease 
(increase) by 1 as kj kj + 1 (kj — > kj — I). If the result- 
ing kj becomes 0, the above process is rejected and take 
the same process again. This process is repeated until 
Kt = Kq is fulfilled. Sort the final set {fci, fe, • • ■ , k^} 
with Kt = Yli=i ki t° satisfy k\ > k 2 > ■ ■ ■ > k^ and 
set ki to be the degree of the node i. Now link the nodes 
2, 3, • • ■ , kx + 1 to the node 1 to make the degree of the 
node 1 be k\. Next, link the nodes k% + 2, fci + 3, • • ■ , 
k% + &2 to the node 2 to make the degree of the node 2 



be &2- Repeat the same processes until no isolated node 
remains. 

Figs. Ola) and Cfb) show snapshots of typical net- 
works generated by the above procedure with 1000 nodes. 
Fig. [Jc) shows the degree distribution of TSFNs con- 
structed with 7 = 2.4 and 4.3. To generate TSFNs we 
use N = 10 5 . The data show a good agreement between 
the exponents 7 obtained from the degree distribution, 
P(k) ~ fc~ 7 , of the constructed network (symbols) and 7 
of the original distribution f(k) ~ fc -7 (solid and dashed 
lines). As shown in Figs, D^a) and^b), the topological 
structure of TSFNs with 7 > 3 is completely different 
from that with 7 < 3. We find relatively small degree 
fluctuation in the networks with 7 > 3 (see Fig. 0(a)). 
The small fluctuation causes almost homogeneous topo- 
logical properties. Thus, we expect that the dynamical 
properties on the networks with 7 > 3 follow the mean- 
field expectations, i.e., the non-zero survival probabilities 
of a lamb in N — > 00 limit. However, for network with 
7 < 3, there exist many short branches and a big hub 
having most links. The dynamical properties of a given 
system on this kind of network can be dominated by this 
hub. Therefore, as we shall show later, the lamb can eas- 
ily hide from the lions by moving to a different branch 
through the hub. 

Now, let's discuss the dynamical properties of the diffu- 
sive capture process on the constructed TSFNs. Initially 
the lamb and n lions are randomly distributed to the 
nodes on the networks. In the simulations, we consider 
two updating rules to determine the position of the ran- 
dom walking lamb and lions: 1) parallel update and 2) 
random update. In the parallel update rule, we move all 
lions and a lamb at the same time, then check if the lamb 
meets the lion(s). In this algorithm, the lamb and the li- 
ons should be initially separated by 2m (m = 1, 2, • • • ) 
steps. Otherwise, the lamb never meets the lion mov- 
ing on TSFNs. In the random update rule, we choose 
one among the lamb or lions at random. The randomly 
chosen lamb or lion takes a random walk and then check 
if the lamb meets the lion. The unit time is defined by 
(n + 1) such random movements. We find that there is 
no significant difference between dynamical behaviors for 
each update rule. We also find that the dynamical be- 
haviors are not varied critically if n is varied. Therefore, 
we will mainly show the simulation results for 1 lamb 
and 10 lions in this article. All the numerical results are 
obtained by averaging over 100 network realizations and 
1000 histories for each network realization. 

Fig. El shows the plot of the average life time (T) on 
TSFNs with various 7 against N using the random up- 
date rule. (T) increases with N for any 7. The lines 
represent the best fits of the relation (T) ~ N a to the 
data. From the fits, we find that (T) follows the power- 
law (T) - N a with a = 1 for 7 > 3, and a = 0.7 for 
7 = 2.4 (dotted line). But, as shown in the Fig. EJ (T) 
for 7 = 2.4 deviates substantially from the power-law re- 
lation (T) ~ N - 7 as increasing N (square). The inset 
shows the results obtained by using the parallel update 
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rule. The data obtained from the parallel update rule 
show nearly the same dynamical behaviors as those ob- 
tained from the random update rule. 

We also measure the survival probability on TSFNs 
with various 7 and N. Fig. [31 shows S(N,t) for 7 = 4.3 
(a), 3.0 (b) and 2.4 (c) from the random update rule. As 
is shown in Fig. [3J the survival probability for 7 > 3 
exponentially decays as S(N,t) ~ exp(— t/r(N)), how- 
ever, S(N, t) for 7 < 3 manifests a substantial devi- 
ation from an exponential decay. If we assume that 
S(N,t) = S e- l l\ then we obtain (T) ~ t(N) from 
the relation (T) = J™ t(-dS(N,t)/dt)dt. Therefore, the 
characteristic time r of S(N,t) scales in the same man- 
ner as (T), and the scaling t/N collapses S(N,t)'s for 
various N into one curve when 7 > 3 [see Fig. O^a) , (b)] . 
The relation t(N) ~ N implies that S(N,t) becomes a 
nonzero constant in the limit N — > 00. Therefore, the 
exponentially decaying S(N, t) for 7 > 3 can be regarded 
as the finite size effect. However, S(N,t) for TSFN with 
7 = 2.4 shows the long-living tail [Fig. which de- 

creases more slowly than exponential decay. Hence, we 
expect 



lim lim 5(7V,t) = { l0ng - livingtail S 7<3 ; (2) 
t^oow^oo v ' | finite (7 > 3) . w 



The origin of 7 dependent behaviors of S(N,t) will be 
discussed in Sec. IV. 

We now discuss the origin of the long-living tail of 
S(N,t) for 7 = 2.4. One possible explanation of the 
slowly decaying tail of S(N,t) can be understood from 
the topological properties of TSFNs. As shown in Fig. 
(b), the biggest hub in the TSFN with 7 < 3 has many 
short branches. Hence, if the lamb hides from the lions by 
moving into a different branch, then the lion should visit 
almost all branches to search for a hidden lamb. Thus, 
the hidden lamb can live longer than the one in the same 
branch with the lion. Therefore, S(N,t) on TSFN with 
7 < 3 can have the long-living tail. In order to verify 
this argument, we count the number of capture events 
N(L, T) with the life time T among 10 6 different histo- 
ries on TSFN with 7 = 2.4 (Fig. 0J) by using the random 
update rule. At the same time we also measure the max- 
imum separation L between the lamb and the lion during 
random walks. Initially we set the shortest path length 
I between the lamb and lion to be 2. We use the TSFN 
whose average diameter d is d ~ 4. If L ~ d, then the 
lamb and the lion have once been in different branches 
during random walks. The data in Fig. ^explicitly shows 
that N(L, t) for L = 2 decays fast (the maximum life 
time of a lamb is T max ~ 10 2 ). In contrast, N(L,t) for 
L = 4 shows the long-living tail up to T max ~ 10 4 . From 
the relation S(t) oc — J J2 L N(L, t)dt, one can easily see 
that the long-living tail of S(t) is caused by the diffusive 
capture process with L = 4. 



III. DIFFUSIVE CAPTURE PROCESS ON 
LSFNS 

The dynamical behavior of one lamb and one lion 
model with the random update rule on LSFNs has al- 
ready been reported in Ref. 29j . We measure the S(N, t) 
and (T) for various n's with the random and the parallel 
update algorithms. We check that (T) and S(N,t) scale 
in the same manner with those reported in Ref. |29|| . 
Therefore, we will report our results briefly in this sec- 
tion. To generate the LSFNs of sizes N = 10 3 ~ 10 6 , 
we use the static network model suggested by Ref. j3]| 
with (k) = 4. In Fig. we display S(N, t) for n = 10 
on LSFNs with various 7(=4.3, 3.0, and 2.4). S(N, t) ex- 
ponentially decays as S(N,t) ~ exp(— t/r) for all 7's. 
As discussed in Sec. II, by assuming that S(N, t) ~ 
exp(— t/r), we obtain the relation (T) ~ r(N). Using 
the value of a's obtained from the best fit of (T) ~ N a 
to the data (which are not shown), we scale the time as 
t/N for 7 > 3 and as t/N - 55 for 7 = 2.4. As shown in 
Fig. Efa), S(N,t) for 7 > 3 collapses well into a single 
curve. This results implies that t(N) ~ N for 7 > 3. 
However, as we can see in Fig. EJc), S(N,t) for 7 = 2.4 
does not collapse into a single curve by the simple scal- 
ing. Note that 7 = 3 is marginal below which (k 2 ) di- 
verges. Thus, the logarithmic correction for the scaling 
of t can be expected for 7 = 3. In Fig. |SJb), we scale the 
time as t/N(logN)- ls and find that S(N,t) for 7 = 3 
collapses into one curve. This means that r scales as 
r ~ iV(log V) -1 - 6 . Therefore, we expect 



lim lim S(N, t) = < )' { 

t^ooA^oo 1 finite (7 > 3) , 



(3) 



where Too is the asymptotic value of t, and depends on 
7. Eqs. (0) and show one of the evidences that the 
dynamical behavior of the diffusive capture process is 
related to (fc 2 ) 



IV. RELATION BETWEEN DEGREES AND 
CAPTURE EVENTS 

In the previous sections, we show that the underlying 
topology significantly affects the dynamical behaviors of 
S(N, t) and (T) of a lamb. For the quantitative descrip- 
tion of the relationship between the structure of the net- 
works and the dynamical properties of lamb-lion model, 
we count the capture events occurring at the nodes with 
degree k. Fig. El^a) shows the normalized number of 
capture events n(k) on LSFNs with 1-lamb and 10-lions. 
As we can see, n(fc)'s both on LSFNs and TSFNs with 
7 = 4.3 decay as n(k) ~ k~ a and the obtained a is 
(7 ~ 2.3. However, n(k) measured on LSFNs with 7 = 2.4 
decreases with a ~ 0.4 for k < 10 3 and increases with 
a ~ -1.35 for k > 10 3 . On TSFNs with 7 = 2.4, we 
also find a ~ 0.4 for 10 2 < k < 10 3 and a ~ —1.9 for 
k > 10 3 (see Fig. Elb)). The decreasing behavior of n{k) 
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for moderate k can be explained by the following anal- 
ysis. In the steady-state, the probability that a random 
walker visits a node i with degree ki is given by |24j | 

P°°=Ak u (4) 

where A = 1/ fcj. Therefore, the normalized capture 
events simply can be written as 

N 

»(*) = E Pi^kPrki = A 2 k 2 NP{k). (5) 

»,j=l 

If P(fc) ~ fc~ 7 , then we obtain 

cr = 7-2. (6) 

Since we do not assume any topological property of un- 
derlying networks, Eq. © is valid for any SFNs and the 
obtained values of a from the simulations for k < 10 3 is 
consistent with Eq. ©. Note that the n(k)'s anoma- 
lously increase for k > 10 3 both on the LSFNs and 
the TSFNs with 7 = 2.4. This anomalous behavior of 
n(k) can be explained by the rarity of the dominant 
hubs in complex networks, i.e., the P(k) for k ps k max 
does not change much as k varies. Here k max is the 
degree of the biggest hub. Thus, if we assume that 
P{k Ri k max ) ~ const. , then n(k) increases as n(k) ~ k 2 
or 

a = -2, (7) 

for k m k max . The obtained a(k m fc maa; )'s for 7 = 2.4 
from the numerical simulations are —1.35 on LSFN and 
— 1.9 on TSFN. The value of a(k w k rnax ) obtained 
from the simulation on LSFN with 7 = 2.4 is consid- 
erably smaller than Eq. (JJJ. This disagreement might 
be originated from our assumption that P{k rj k max ) 
does not depend on k. For example, if we assume that 
P{k rs k max ) ~ fc~ e with e < 2, then tr becomes 
cr = e - 2(< 0). Eqs. JSJ) and 0, thus, provide the 
topological origin of the gathering behavior of random 
walkers at the dominant hubs. It implies that the lamb 
and the lion have a strong tendency to move into big 
hubs, i.e. the hubs collects the lamb and the lion. There- 
fore, the hubs in the networks with strong heterogeneity 
become the effective attractors of random walking lamb 
and lions. 

Another important consequence of Eq. (J5J is that the 
dynamical properties of lamb-lion model depend on (fc 2 ) . 
Note that from Eq. (JSJ the total number of capture event 
Ntot can be expressed as 

AW cx J>(fc) ~ J> 2 P(fc) ~ (fc 2 ) . (8) 

k k 

If 7 < 3 then N to t also diverges. This means that the 
lions on the networks can capture the lamb no matter 



how the lamb is trying to escape from the lions. There- 
fore, S(N,t) -» in the limit A -> 00. For 7 > 3, N tot 
converges to a finite value and the lions can catch the 
lamb with a probability p < 1, thus S(N, t) ^ for large 
A. The expectations agree well with our simulation re- 
sults shown in Sec. II and III. Based on Eq. JHJ, we also 
expect that the long-living tail of S(N, t) on TSFN with 
7 < 3 eventually goes to 0. 



V. SUMMARY AND DISCUSSION 

We investigate the diffusive capture process on two dif- 
ferent kinds of SFNs. The diffusive capture process can 
be mapped into the diffusing lamb and lion model. We 
measure S(N, t) and (T) of a lamb. We find that on both 
TSFNs and LSFNs S(N,t) of a lamb decays exponen- 
tially for 7 > 3. Simple scaling arguments predict that 
the exponentially decaying S (A, t) for 7 > 3 is originated 
from the finite size effect. This result recovers the mean- 
field expectation S(N,t) —* finite in the limit N — > 00. 
However, we find that S(N, t) on TSFNs with 7 < 3 has 
the long-living tail. By measuring the number of capture 
events with life time T we find that the long-living tail 
is caused by enormously many short branches in TSFNs. 
On LSFN with 7 < 3, we find that S(N, t) decays ex- 
ponentially and no lamb can survive in thermodynamic 
limit. 

In order to understand the non-mean-field behavior of 
S(N,t) for 7 < 3, we study the normalized number of 
capture events n(k) both numerically and analytically. 
From the numerical simulations, we find that n{k) de- 
creases as a power-law, n(k) ~ k~ a . The simple ana- 
lytic derivation of n(k) provides the relation a = 7 — 2, 
which shows a good agreement with numerical results for 
moderate k. However, we find that n(k) anomalously in- 
creases for large k in the simulations. This means that 
the most of the lambs are captured by lions at the hubs. 
Hence, those hubs become the effective attractors. By 
assuming the rarity of dominant hubs and P(k) w const. 
for large k(~ k max ), we show that the n(k) can increase 
for large k. From the analytic expectation of n(k), we 
also show that Ntot ~ (^ 2 )- This implies that if the 
fluctuation in degree is large enough (7 < 3) then the 
dynamical behaviors of the diffusive capture process de- 
viate from the mean-field expectations by emerging ef- 
fective attractors of random walkers. From the analytic 
expectation of n{k) and A tot , we also expect that the 
long-living tail of S(N, t) observed on TSFN with 7 < 3 
(see Eq. (J5J) will eventually go to due to the diverging 
(k 2 ). 

This work is supported by Korea Research Foundation 
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FIG. 1: Schematic structures of TSFNs with 7 = 4.3 (a), 
and 2.4 (b) constructed by the method introduced in Sec. II. 
The plot of the degree distribution of the TSFNs with 7 = 2.4, 
and 3.0 (c). Network size is N = 10°. The degree distribution 
P(k) ~ fc~ 7 of the constructed networks follows the original 
distribution f(k) ~ fc~ 7 well. 




FIG. 2: Average life time (T) on TSFNs with 7 = 4.3, 3.0 
and 2.4. The lines represent the best fits of (T) ~ A*" 1 to data. 
The obtained exponents are a ~ 1.1 (solid line) for 7 = 4.3, 
a — 1.0 (dashed line) for 7 = 3.0 and a = 0.7 (dotted line) 
for 7 = 2.4. For 7 = 2.4 (T) substantially deviates from 
(T) ~ iV ' 7 line when N = 10 6 . 
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FIG. 3: The scaled plots of the survival probability 5 , (A r , t) 
on TSFNs with 7 = 4.3 (a), 3.0 (b) and the plot of S(N,t) 
for 7 = 2.4 (c). For the scaling plots, a = 1 which is obtained 
from the fits for 7 > 3 in Fig. |5]is used. 



8 




FIG. 4: Plot of N(L, T) on TSFN with 7 = 2.4 and N = 10 4 . 
The initial separation I of the lamb and the lion is 2. L is the 
maximum separation between the lamb and the lion during 
random walks (symbols). 
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FIG. 5: The scaled plots of the survival probability S(N, t) 
on LSFNs with 7 = 4.3 (a), 3.0 (b), and 2.4 (c). To scale the 
time t, a's are obtained from the fits of (T) ~ N a to the data 
(which are not shown in this article). 
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FIG. 6: Plot of capture event distribution n(k) on LSFN 
(a) and TSFN (b). The lines stand for the relation n(k) ~ 
k~ a . The solid lines correspond to a = 0.4 and the dotted 
lines represent a = 2.3. The data for k < 10 3 (a) and for 
10 2 < k < 10 3 (b) show good agreements with the relation 
a — 7 — 2. The dashed lines correspond to a = — 1.35 (a) and 
a = —1.9 (b) representing the anomalous increases of n(k)'s 
for large k. 



